Hidden Caldeira-Leggett dissipation in a Bose-Fermi Kondo model 
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We show that the Bose-Fermi Kondo model (BFKM), which may find applicabihty both to certain 
dissipative mesoscopic qubit devices and to heavy fermion systems described by the Kondo lattice 
model, can be mapped exactly onto the Caldeira-Leggett model. This mapping requires an ohmic 
bosonic bath and an Ising-type coupling between the latter and the impurity spin. This allows us 
to conclude unambiguously that there is an emergent Kosterlitz-Thouless quantum phase transition 
in the BFKM with an ohmic bosonic bath. By applying a bosonic numerical renormalization group 
approach, we thoroughly probe physical quantities close to the quantum phase transition. 
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The Bose-Fermi Kondo model (BFKM) (or equiva- 
lently the spin-boson- fermion model), originally intro- 
duced by Si and coworkers^^ and by Sengupta^ to de- 
scribe peculiar quantum critical behaviors in heavy- 
fermion Kondo lattice systems^, involves a single impu- 
rity spin being coupled both to a bosonic bath and to a 
fermionic bath. Resulting from the nontrivial competi- 
tion between these distinct baths, a rich phase diagram 
is known to emerge from this model (See, e.g., Refs. 
and0). More generally, a great interest is currently de- 
voted to the understanding of the Kondo entanglement 
breakdown mechanism due to the presence of extra (here, 
bosonic) quantum fluctuations resulting in striking quan- 
tum phase transitions^. In this Letter, we revisit the case 
where the impurity spin S* = 1/2 is coupled to an ohmic 
bosonic bath with a continuum spectrum - ohmic means 
that the bosonic correlation function in time t decays as 
- through an Ising coupling. The anisotropic Hamil- 
tonian under consideration thus takes the form: 

H = hS, + H,f +Hsb, (1) 
^^./ = %(*l(0)*T(0)^++h.c.) 

/oc 
dx ^l{x)id,^„{x), (2) 

y/2Kb 47r J^^ 

where h \s a. magnetic field, '^cr{x) and <& represent the 
fermionic and bosonic fields; Vh is the velocity of the 
bosons and K^"^ ^ the typical coupling between the 
bosons and the impurity spin; Recall that K^^ ~ 
means no coupling between the impurity spin and the 
bosonic environment. In Eq. (01, Vf is the Fermi velocity. 
Without loss of generality, the one-dimensional character 
of the Hamiltonian Hsj can be viewed as a result of the 
point-like character of the impurity and the rotational 
symmetry. Besides, we omit the Ising part of the Kondo 
coupling Jz due to its minor effect (see the last page). 

Our interest in the model is also motivated by the 
fact that this model may be realized in mesoscopic dis- 
sipative setups involving qubits, as pointed out by one 



of us recently''. More precisely, the impurity spin can 
embody the two allowed charge states of a big metallic 
grain close to a given degeneracy point and h being pro- 
portional to the gate voltage measures deviations from 
this degeneracy point^. The conduction electrons stand 
for the electrons both in the metallic grain {'^i) and in 
a nearby reservoir electrode (^'-f), and the J± (Kondo) 
term denotes the tunneling process of an electron from 
lead to grain that flips (through the raising operator S+) 
the charge state of the grain, and vice-versa. Note that 
here the spin index cr = ("f , i) is completely artificial and 
refers to the position of an electron in the structure (lead 
or grain) . The original spin of the electrons are assumed 
to be polarized ("spinless electrons") due to the appli- 
cation of a strong magnetic field. The bosons represent 
the electromagnetic noise in the gate voltage stemming 
from the finite resistance R in the gate lead, and the spin- 
boson coupling reflects the effect of the voltage noise on 
the charge fluctuations of the grain. Here, Kb = Rk/2R 
with Rk = 27r?i/e^ the quantum of resistance^i. Two of 
us have extended this model to one-dimensional reservoir 
leads being embodied by a Luttinger liquid behavios^ and 
also to two strongly capacitively-coupled large quantum 
dots (which can be reduced to a single charge qubit) 

As a first step, the model of Eqs. (|1I3|I has been exten- 
sively studicd^'^'^ by using a perturbative renormaliza- 
tion group (RG) approach. A quantum phase transition 
is discovered at the critical value of i^gLii: 

(Kb)-' = J^/nvf = 2AK, (4) 

where ujc is the high energy cutoff in the theory and A — 
J±u}c/'2nvf. The RG fiow equations suggest that the 
phase transition is of the Kosterlitz-Thouless (KT) type. 
On the other hand, a more rigorous approach is required 
in order to unambiguously prove the KT transition as 
well as to scrutinize the evolution of physical quantities 
in the vicinity of the quantum phase transition. In this 
Letter, we explore an exact mapping of the BFKM in 
Eqs. (|1I3|I : We properly demonstrate that this model can 
be mapped onto the Caldeira-Leggett (CL) modelM with 
the effective dissipation strength 

a = l + iAKb)-\ (5) 
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and therefore the two models belong to the same class 
of universality. It immediately follows that there is a 
KT quantum phase transition in the BFKM, separating 
a Kondo phase and an unscreened spin phase. Moreover, 
it irrefutably demonstrates the important Eq. (4). Even 
though the Bethe-Ansatz methodic can be exploited in 
the Kondo region, it breaks down in the vicinity of the 
quantum phase transition^^ (which must be identified 
as the antiferromagnetic-ferromagnetic transition in the 
anisotropic Kondo model) . We resort to a bosonic numer- 
ical RG (NRG) technique applying to the CL model^- . Of 
interest to us here is the spin magnetization at tempera- 
ture T = 0, (Sz), as well as the (local) spin susceptibility 
Xioc{T) versus T. In the mesoscopic realizations, (Sz) = 
{Q) — 1/2, where Q is the charge operator on the (large) 
dot, can be measured with very high precision^'^. 

CL mapping. — We first bosonize the fermions in 
Eq. Q as: *|(x)5'|(a;) = [wc/(27rw/)] exp(z\/2(p(a;)) and 
«'|(a:)^'|(a;) - ^[{x)^ ^{x) = dMx)/iV2n); A and Uc 
are defined in Eq. (4). By dropping the trivial contribu- 
tion in the charge sector, the Hamiltonian H^f becomes 



h.c. 



/oo 
dx\dxip{a 



,(6) 



and the Hamiltonian |^ becomes that of an impurity 
spin coupled to two bosonic baths. Now we show that 
in the case of an ohmic bosonic bath in Eq. (3) one 
can envision to combine the two bosonic baths and even- 
tually map the model onto a well-known CL model. 
The local action for the boson fields (p{x) reads S*^""^ ~ 

kn|</'o(w„)v3o(-^n), where ip^ = ip(x = 0) and 
LiJn is the bosonic Matsubara frequency. Such a form 
implies that ip{x) is of the ohmic type. If the bosonic 
bath $(a;) in Eq. (|3Jl is also ohmic and thus embodied 
by the local action S^^" = ^J2u„ |cj„|$o(wn)$o(~t^«) 
with $0 = = 0), we can then employ linear com- 
binations of (fio and $0 to simplify the problem. (Note 
that starting from a non-ohmic bosonic environment^® 
instead, one would replace |w„| in by \ujn\^ with 
s ^ 1; the two bosonic baths thus would have different 
dynamics which would hinder the straightforward combi- 
nation.) By performing the unitary transformational — 
exp{-i^oSz/V^Kb} to Eqs. Ql, 01 and ©, we can ab- 
sorb the spin-boson coupling term in Eq. Q into the 
transverse Kondo term A in Eq. It is then convenient 
to introduce the symmetric and antisymmetric bosonic 
combinations ips = (2a)"^/^[V2(y3o + {2Kb)~^^^ ^0] and 
(pa = (2q!)- [(27^^) -1/2(^0 - V2^o], where the crucial 
parameter a has been defined in Eq. |(SJ). The impurity 
spin gets only coupled to the ips mode and it is enough 
to write down the (ySs-action 



\uJn\(Ps{l^n)fs{-l^n) 
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2^ 



Conceptually, we can visualize S*,^^ as the action linked to 
an Hamiltonian H^p^ with the effective velocity = uJcCi 
(a is the short-distance cutoff which will be used be- 
low) and modeling a single bosonic bath with the dis- 
sipative parameter a coupled to the impurity spin. The 
link with the CL model of a two-level system with ohmic 
dissipationi^ becomes clear when applying the unitary 
transformation = exp{i\/2aips{0)Sz} to turn the A 
term into a transverse magnetic field ASx, leading to 

H = Hs- V2^v,d^ip,{0)Sz + / dx[dx^s{x)f,{8) 

where Hs — hSz + and a is given in Eq. (|SJ). 

This mapping is crucial because it irrefutably reveals 
an originally hidden correspondence between the BFKM 
of Eqs. (1-3) and the CL model in which the breakdown 
of the Kondo physics (quantum coherence) due to dis- 
sipation is clearly established. For example, a KT type 
transition was predicted in Ref. 12 in the context of ohmic 
dissipation and this has been solidly confirmed through 
the NRG, either by refermionizing the bosonic degrees 
of freedom and then mapping the CL model onto the 
anisotropic Kondo model^^, or directly by resorting to 
the bosonic NRCU^. For A — > the KT transition is 
known to occur at^^ a — ac = 1 implying in the BFKM 
a critical (iff,)"^ from Eq. jSJ: since A 0, the 
Kondo energy scale Tx characterizing the emergence of 
a bound state between the local moment and the con- 
duction electrons vanishes and thus the local moment 
remains unscreened whatever the coupling between the 
dissipative mode dx^{0) and the local spin. 

To understand more deeply the connection between 
the CL model and the Kondo phenomenon, we can pro- 
ceed from Eq. (8) along the lines of Ref. and pre- 
cisely recover the anisotropic Kondo model. The effective 
Kondo parameters here are given by J± = nVgA/ujc and 
Jz — 4:iTVs{l — y/a). Thus, one can clearly distinguish two 
phases; a Kondo (delocalized) realm when Jz > —\J±\ 
and a ferromagnetic Kondo phase which embodies an 
unscreened (localized) moment when Jz < — |Jj^|. This 
enables us to rigorously predict a KT phase transition 
in the BFKM when Jz — — |J_l| which implies a crit- 
ical Qfc = 1 + {Kb)c^ /4: = 1 + A/(2wc) when assuming 
< A ^ LUc', this reinforces the intuition gained from the 
perturbative RG analysis of the BFKM''. Since Bethe- 
Ansatz calculationsi^ can only be applied to this model 
when a < we seek to apply the (bosonic) NRG on 

Eq. (8) to probe (Sz) versus h and xioc{T) very close to 
Uc- Those quantities have not been studied in Ref. [l^ 

NRG endeavors close to the transition. — The strength 
of the NRG lies in its nonperturbative nature and the 
ability to resolve arbitrarily small energiesi2i22i2i. This 
allows to provide important information in the phase 
transition region. For convenience, the Hamiltonian 
can be rewritten as^S 



(7) 



H = Hs + ojc 



de 



+ e ala 



(9) 
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N h 
FIG. 1: (color online) Left: Flow of (Sz) with the iteration 
N. The dashed line is for h = —0.0001 (unscreened realm). 
All the solid lines correspond to a = 0.6 (Kondo regime), and 
from top to bottom, to h = —10""', —5 x 10""'', — 10~^, —5 x 
10"^-10"^ Right: (S,) as a function of h. The open circles 
sitting on top of the a = 0.6 curve are the Bethe-Ansatz 
resultai^. Parameters are A = 0.01 and ujc = 1. The NRG 
parameters are A = 2, A^s = 100, Af, = 8, and Am = 500. 



Then, we introduce a logarithmic discretization of the 
energy interval [0, 1], i.e., a decreasing set of frequencies 
(1,A"\---,A"",---) with A > 1 and inte ger n. After 
transforming to a semi-infinite chain, H reads^^^ 



Hn 



lim 



N- 



N, 



A 



N 



Hs + ^f^uJc (bo + bl)S, 



N 



n=0 



N~l 

E 

n=0 



On+i -t- n.c.) 



(10) 



(11) 



In Eq. H1U|I the spin gets only coupled to the first (0th) 
site of the bosonic chain, and the remaining part of the 
chain is characterized by on-site energies e„ and hopping 
parameters t„. They satisfy a set of recursion relations 
which can be solved numerically. We refer the reader to 
Ref. for the details of these recursion relations and 
their derivations as well as for the precise definition of 
the bn boson operators. It is important to note that 
both e„ and tn decay exponentially as A~" as a result 
of the logarithmic discretization. This allows us to solve 
the model in Eq. (|10|) in an iterative way: first we di- 
agonalize ifjv exactly, keep the Ns lowest energy levels 
and then use the recursion relation between Hj^^i and 
Hn derived from Eq. (fTT|l to diagonalize Hm+i. Typi- 
cally, the result converges at large iV ~ 25. In contrast to 
the fermionic case, each boson site allows for an infinite 
number of bosons and hence a truncation of the {Nj, -\- 1) 
boson states is necessarjiS^. Throughout this paper, we 
use Ns = 100 and Nf, = 8 (except the 0th site for which 
Nto = 500); those values are approximately those used 
in Ref. 0. We have also used another set of parameters 
{Ns = 150 and Nt = 12) and found a good agreement. 
The quantities {S,){h) and xioc(r) = d{S,){T)/dh\h^o 
are also calculated iteratively. 

NRG Results. — In Fig. [T] (left) we show the flow of the 
ground state expectation value (Sz) as a function of the 
iteration variable N. At small A/wc, for a = 1.4 and for 
any h < 0, after a few iterations we find that (Sz) flows 




FIG. 2: (color online) The local susceptibility xioc(r) as a 
function of T. The open circles are the numerical data for 
T = Tn ~ 1.1A~^ and lines are guides to the eyes. Inset: 
Logarithm of the Kondo energy scale Tk versus a for a close 
to Qc. The solid circles are the numerical data and the dashed 
line represents the fitting curve InT/f = —4.67/(1.15 — a). The 
NRG parameters, A and ujc, are the same as in Fig. 



rapidly to 1/2, which is in agreement with the localized 
(dissipative) phase of the CL modeli^ as well as with 
the perturbative RG analysis of the BFKM"^. This corre- 
sponds to an unscreened local spin or a particle which is 
localized in one of the two allowed levels. A strong dis- 
sipative Ising coupling with the bosonic bath definitely 
hinders the Kondo effect to develop. In contrast, for 
a = 0.6, we show that (Sz) flows to saturating values 
which deviate from 1/2 and fall drastically with the de- 
creasing h. Moreover, (Sz) yields a linear dependence 
versus h around h = which indicates a fully screened 
moment: this is a clear signature of the Kondo effect. 
We have fitted the NRG curve for a = 0.6 (which is suf- 
ficiently far from ac) with Bethe-Ansatz results applying 
to the Kondo reahr.^'^-^^. The dependence of {Sz) on h at 
various a for A — 0.01 and ujc = 1 is displayed in Fig. ^ 
We have carefully followed the progressive destruction 
of the delocalized Kondo regime when approaching the 
quantum phase transition at different applied magnetic 
fields. More specifically, we can infer the Kondo tempera- 
ture Tk from {Sz) versus h based on the following scaling 
argument. The (Kondo) Fermi liquid type state emerging 
for a < Qfc is embodied by a constant (local) spin suscep- 
tibility xioc = d{Sz)/dh which must be identified as2ii^ 
Xioc ~ ^/Tk- For a not too far from ac, we can then ex- 
tract Tk as the width of h for {Sz) changing from -0.45 
to 0.45. This NRG procedure is complementary to that 
from the energy spectrum flowi^. We check that the re- 
sulting Tk obeys InT^: cx l/{ac — Oi) close to ac as clearly 
shown in the inset of Fig. |2]for A = 0.01 and Wc = 1~- 
The fact that Tk decreases exponentially fast close to ac 
is definitely reminiscent of the KT transition. The scaling 
in the inset of Fig. 2 also reveals that for the used NRG 
parameter A = 2, ac = 1.15. In fact, the NRG result is 
only exact in the limit of A ^ 1. Hence, the precise ac 
can be obtained by getting ac for a few A parameters and 
by extrapolating the result to A ^ 1. For A/lJc = 0.01, 
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we reach ac(A — » 1) = 0.98 (~ 1). The effect of the dissi- 
pative coupUng becomes weaker when increasing A. For 
each A, we have repeated the above-mentioned extrapo- 
lation procedure and precisely extracted Uc which basi- 
cally corresponds to the critical bosonic coupling at which 
the Kondo scale strictly vanishes. The resulting ac as a 
function of A obeys ac(A) = 0.974 + 0.537 A/ujc which 
clearly agrees with our expectation from the CL map- 
ping. The formation of a Fermi liquid ground state in the 
Kondo realm is also clearly shown in Fig. |21 through the 
X\oc{T) versus T plot: For a < — 1.15, XiodT) satu- 
rates to a constant for T < Tk whereas for T above Tk 
the Curie's law xioc oc 1/T is well satisfied. At the tran- 
sition point, we clearly observe that the xioc(2^) c>c 1/T 
behavior persists down to the lowest temperatures. 

Finally, we would like to briefly comment on the effect 
of the Ising coupling term Jz between the spin and the 
fermions, originally present in the heavy fermion systems 
but absent in the mesoscopic realizations under consid- 
eration, close to the phase transition point. For a — > 1, 

~ (/3s and thus the second term in the right hand side 
of Eq. © becomes \/2(^J^ - ^Vs)da,ips{'i))Sz. This 
immediately leads to a critical Uc — I + C(A) + 0{Jz)^ 
which is qualitatively consistent with Ref. |25|. 



Conclusion. — We have demonstrated that the BFKM, 
with an ohmic bosonic bath and an Ising-type coupling 
between the latter and the impurity spin, and the CL 
model belong to the same class of universality. This un- 
ambiguously proves the existence of a KT phase transi- 
tion at zero temperature in the former model when tuning 
the Ising dissipative coupling. We have chosen to resort 
to the bosonic NRG to investigate in detail physical quan- 
tities like the spin magnetization {Sz) versus the mag- 
netic field h close to the quantum phase transition where 
Bethe-Ansatz calculations break down. The latter clearly 
reflects two phases: the "localized" phase typical of an 
unscreened moment and the "delocalized" Kondo realm 
characterized by a Kondo energy scale which goes to zero 
exponentially fast close to the KT transition. Compar- 
ing with the recent fermionic NRG of the BFKM^^, the 
bosonic NRG approach can be extended to the case of 
a sub-ohmic bat h^^i^^ which might be relevant for the 
quest of quantum criticality in heavy-fermion systems ». 
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